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C^; ABSTRACT 

o ■ 

a: 

We present the Observations of Redshift Evolution in Large Scale Environ- 
ments (ORELSE) survey, a systematic search for structure on scales greater than 

CN ■ 10 h^Q Mpc around 20 well-known clusters at redshifts of 0.6 < z < 1.3. The goal 

^ ■ of the survey is to examine a statistical sample of dynamically active clusters and 

large scale structures in order to quantify galaxy properties over the full range 
■ of local and global environments. We describe the survey design, the cluster 

Q\ ■ sample, and our extensive observational data covering at least 25' around each 

target cluster. We use adaptively-smoothed red galaxy density maps from our 
wide-field optical imaging to identify candidate groups/clusters and intermediate- 
density large scale filaments/walls in each cluster field. Because photometric 

^ . techniques (such as photometric redshifts, statistical overdensities, and richness 

estimates) can be highly uncertain, the crucial component of this survey is the 
unprecedented amount of spectroscopic coverage. We are using the wide-field, 
multi-object spectroscopic capabilities of the DEep Multi-Object Imaging Spec- 
trograph to obtain 100-200+ confirmed cluster members in each field. Our survey 
has already discovered the CI 1604 supercluster at z ~ 0.9, a structure which con- 
tains at least eight groups and clusters and spans 13 Mpc x 100 Mpc. Here, we 
present the results on the large scale environments of two additional clusters, 
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Cl 0023+0423 at z = 0.84 and RX J1821. 6+6827 at z = 0.82, which highlight 
the diversity of global properties at these redshifts. The optically-selected CI 
0023+0423 is a four-way group-group merger with constituent groups having 
measured velocity dispersions between 206-479 km s _1 . The galaxy population 
is dominated by blue, star-forming galaxies, with 80% of the confirmed members 
showing [Oil] emission. The strength of the H<5 line in a composite spectrum of 
138 members indicates a substantial contribution from recent starbursts to the 
overall galaxy population. In contrast, the X-ray-selected RX J1821. 6+6827 is a 
largely-isolated, massive cluster with a measured velocity dispersion of 926 ± 77 
km s -1 . The cluster exhibits a well defined red sequence with a large quiescent 
galaxy population. The results from these two targets, along with preliminary 
findings on other ORELSE clusters, suggest that optical selection may be more 
effective than X-ray surveys at detecting less-evolved, dynamically-active systems 
at these redshifts. 

Subject headings: catalogues - surveys - galaxies: clusters: general - large-scale 
structure of the Universe 



1. Introduction 

Clusters of galaxies are accurate tracers of the large scale structure in the local Universe. 
Redshift surveys at z < 0.1 (e.g., Geller & Huchra 1989; Einasto et al. 1997; Colless et al. 
2001) and numerical simulations (e.g., Colberg et al. 2000; Evrard et al. 2002; Dolag et 
al. 2006) reveal the filamentary structure of the Universe, stretching between clusters and 
superclusters of galaxies. Clusters form at the nodes of these filaments, growing through 
the continuous accretion of individual galaxies and groups from the surrounding field (e.g., 
Frenk et al 1996; Eke et al. 1998). The precursors of the filaments should be present around 
distant clusters, containing many of the galaxies which will eventually infall into the virialized 
core and form the cluster population observed today. Since galaxy environment should 
change dramatically during the course of vigorous cluster assembly, the large scale structure 
present around high-redshift clusters offers us the unique opportunity to probe, over the full 
range of local environments, the physical effects on galaxies as they assemble into denser 
regions. As a result, we are undertaking the Observations of Redshift Evolution in Large 
Scale Environments (ORELSE) survey, a systematic search for structure on scales greater 
than 10 h^Q Mpc around 20 known clusters at z > 0.6. The survey covers 5 square degrees, 
all targeted at overdense cluster regions, making it complementary and comparable to field 
surveys such as DEEP2 (Davis et al. 2003) and COSMOS (Scoville et al. 2007; Koekemoer 
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et al. 2007). 

While superclusters have been studied at low redshift (Davis et al. 1980; Postman, Geller 
& Huchra 1988; Quintana et al. 1995; Small et al. 1998; Rines et al. 2002; Fadda et al. 2008; 
Porter et al. 2008) and at z m 0.2 - 0.6 (Kaiser et al. 1999; Gray et al. 2002; Kodama et al. 
2001, 2003, 2005; Ebeling et al. 2004; Tanaka et al. 2007a; Kartaltepe et al. 2008), large scale 
structures at high redshift (z > 0.6) are just now being explored (Kodama et al. 2005; Nakata 
et al. 2005; Tanaka et al. 2006, 2007b; Swinbank et al. 2007; Gal et al. 2008; Fassbender 
et al. 2008; Gilbank et al. 2008; Patel et al. 2008). Previous studies at these redshifts have 
largely focused on the central regions (< 2 Mpc) of massive clusters. These studies have 
revealed strong evolution in the cluster galaxy population which includes (1) an increase in 
the fraction of blue, star-forming, late-type galaxies with redshift, implying that early- types 
are forming out of the excess of late- types over the last ~ 7 Gyr (e.g., Butcher & Oemler 
1984; Ellingson et al. 2001; Dressier et al. 1997; van Dokkum et al. 2000, 2001; Lubin et 
al. 2002); (2) a larger fraction of post-starburst ("K+A") galaxies in the cluster versus field 
environment, indicating strong star formation activity in the recent past (e.g., Dressier et al. 
1999, 2004; Balogh et al. 1999; Tran et al. 2003); (3) a deficit of faint, passive, red galaxies, 
suggesting that a large fraction of these galaxies in present-day clusters has moved onto the 
red sequence relatively recently as a result of a truncation of star formation (e.g., Smail et 
al. 1998; van Dokkum & Franx 2001; Kodama et al. 2004; De Lucia et al. 2004, 2007; Tanaka 
et al. 2005, 2007b; Koyama et al. 2007); and (4) an increase in the overdensities of X-ray 
and radio sources with cluster redshift, indicating enhanced starburst /nuclear activity in the 
past (e.g., Best 2003; Cappelluti et al. 2005; Eastman et al. 2007; Kocevski et al. 2008a). 

All of these findings imply significantly increased star-forming, starburst, and nuclear 
emission in the past and raise the question - what physical mechanisms associated with 
the cluster environment are responsible for the suppression of star formation and nuclear 
activity and the transformation of gaseous, disk galaxies into passive spheroids? Several 
mechanisms have been suggested, including galaxy harassment (Moore et al. 1998), ram 
pressure stripping (Gunn & Gott 1972), starvation (Larson et al. 1980), and merging (Mihos 
1995, 1999). Most of these processes are associated not with the densest cluster regions, but 
rather with the infall regions and lower-density environments far from the cluster cores. 

Studies at low-to-moderate redshift suggest that the non-cluster processes play a pivotal 
role in driving galaxy evolution well before the galaxies reach the cluster cores. Data from 
the Sloan Digital Sky Survey (SDSS) and the 2dF Survey indicate a sharp transition between 
galaxies with field-like star formation rates (SFRs) and galaxies with low SFRs comparable 
to that of galaxies within the cluster cores (Lewis et al. 2002; Gomez et al. 2003; Goto et 
al. 2003a,b; Balogh et al. 2004). The density at which this transition occurs (log £ ~ 1 
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Mpc -2 ) corresponds to the cluster virial radius, well outside the core region on which most 
studies have been focused. The fact that low SFRs and passive (gas-less) spiral galaxies 
are observed well beyond the virial radius rule out severe processes, such as ram pressure 
stripping or galaxy-galaxy merging, as being completely responsible for the variations in 
galaxy properties with environment (Lewis et al. 2002; Goto et al. 2003b). Similarly, studies 
on large scales (clustocentric distances of > 4 hjQ Mpc) at z ~ 0.4 indicate that galaxy 
colors change sharply at relatively low densities and that the morphological mix, at a given 
density, is independent of cluster radius, implying that galaxy transformation occurs outside 
the cluster core and is physically associated with the infalling filaments and the chains of 
smaller galaxy groups within them (Kodama et al. 2001; Treu et al. 2003). Consequently, 
these groups serve as a pre-processing phase in the evolution of cluster galaxies (Kodama et 
al. 2001, 2003, 2004; Gray et al. 2002; Treu et al. 2003; Bower & Balogh 2004; Mihos 2004). 

These trends continue to redshifts approaching z ~ 1 where there are indications of 
"downsizing" with galaxy evolution accelerated in high-density and high-mass (cluster) re- 
gions, compared to lower-density and lower-mass (group) regions. Studies of large scale 
structures at z ~ 0.5 — 0.8 suggest that the galaxy population in clusters is more evolved, 
with the faint-end of the red sequence being more developed and the red-sequence galaxies 
showing no signs of current or recent star formation. In contrast, groups and lower-density 
environments have a stronger deficit of the faint red-sequence galaxies and clear signs of 
star- format ion activity (weak [Oil] emission and/or strong B.5 absorption) in the existing 
red-galaxy population (De Lucia et al. 2004, 2007; Kodama et al. 2004; Tanaka et al. 2005, 
2005; Koyama et al. 2007). As such, these results confirm that, at these epochs, much of the 
activity is taking place in lower-density, lower-mass environments. 

This environment-activity level connection also extends to the extremes of active galac- 
tic nuclei (AGN) and starburst galaxies. From studies covering a wide redshift range (0.2 < 
z < 1.2), excesses of 24 /im and X-ray sources are observed preferentially in intermediate- 
density, not high-density, regimes, consisting of groups, cluster infall regions, and filaments 
(Cappelluti et al. 2005; Gilmour et al. 2007; Marcillac et al. 2007, 2008; Kocevski et al. 
2008b). These results, as well as those described above, imply that the environmental pro- 
cesses which induce nuclear activity, truncate star formation, and change galaxy properties 
are not necessarily driven by cluster specific mechanisms, like ram pressure stripping by the 
hot intracluster medium or harassment (i.e., truncation of the galaxy halo) by the cluster 
tidal field. Thus, it is essential to look well beyond the regions of traditional study (the 
cluster cores) to find answers to our questions concerning galaxy evolution. The ORELSE 
survey is designed specifically to target these largely-unexplored regions and answer these 
questions by correlating multi-wavelength (radio, optical, infrared, and X-ray) photomet- 
ric data with galaxy kinematics and spectral features in a statistical sample of large scale 
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structures at redshifts approaching unity. 

The first large scale structure detected in the ORELSE survey was the CI 1604 su- 
percluster at z ~ 0.9 which includes two massive clusters CI 1604+4304 at z — 0.90 and 
CI 1604+4321 at z = 0.92 originally detected by Gunn, Hoessel & Oke (1986) and further 
studied by Oke, Postman & Lubin (1998). Wide-field imaging and spectroscopy as part of 
the ORELSE survey has revealed a more complex, massive structure, containing at least 
eight groups and clusters and spanning 13 Mpc x 100 Mpc (Lubin et al. 2000; Gal & Lubin 
2004; Gal, Lubin & Squires 2005; Gal et al. 2008). Our extensive spectroscopic data on 
the CI 1604 supercluster (over 400 confirmed members) demonstrate that comprehensive 
redshift surveys, like ORELSE, are essential for understanding galaxy and cluster evolution. 
Specifically, (1) the entire structure size in redshift space is equivalent to typical photometric 
redshift errors (Margoniner & Wittman 2008); (2) superpositions of groups/clusters mean 
that mass measures based on weak lensing signal, richness, or X-ray luminosity are highly 
uncertain; (3) cluster velocity dispersions based even on traditionally large numbers of galax- 
ies (i.e. 20-40) can be substantially overestimated due to outliers (Gal et al. 2008); and (4) 
the expected overdensities of radio and X-ray sources are small, making the identification 
and study of individual active galaxies impossible (Kocevski et al. 2008a). 

In this paper, we build on our studies of the CI 1604 supercluster. We present the 
experimental design of the full survey (§2) and the photometric and spectroscopic results 
on two additional target clusters, CI 0023+0423 (hereafter CI 0023) at z = 0.84 and RX 
J1821. 6+6827 (hereafter RX J1821) at z = 0.82, which highlight the significant diversity of 
structure properties at these redshifts (§4 and 5). Throughout the paper we use a cosmology 
with H = 70 /ifo 1 km s" 1 Mpc" 1 , Vt m = 0.3 and fi A = 0.7. 



2. The Survey Design 

The ORELSE survey is a systematic search for structure on scales greater than 10 H-jq 
Mpc around 20 well-studied clusters at redshifts of 0.6 < z < 1.3 (see Tabled]). The goal 
of the ORELSE survey is to study the large scale structure around a statistical sample of 
high-redshift clusters. Numerical simulations indicate that 50% of all M > 3 x 10 14 hj M Q 
clusters at z ~ 0.8 have two or more companion clusters of equal or greater mass within 100 
co-moving Mpc (i.e., a configuration similar to the CI 1604 supercluster), and 70% have had 
its most recent large merger (M f ina i / 'M init iai > 1.2) within the last ~ 1 Gyr (Colberg et al. 
2000; Evrard et al. 2002; Cohn & White 2005). Thus, our survey should identify at least 
ten structures like the CI 1604 supercluster and detect significant formation activity in the 
majority of clusters in our sample. 
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2.1. The Cluster Sample 

The cluster sample is chosen from various surveys (Oke, Postman & Lubin 1998; Rosati 
et al. 1998; Vikhlinin et al. 1998; Gladders & Yee 2000; Gioia et al. 1999, 2003; Stanford 
et al. 2002; Blanton et al. 2003; Hashimoto et al. 2005; Henry et al. 2006; Pierre et al. 
2006; Maughan et al. 2006) with different detection techniques (optical, X-ray, radio) to 
reduce any bias associated with sample selection. Table [1] lists the cluster coordinates, 
redshift, detection technique, availability of high-quality (Chandra and/or XMM) X-ray 
data, status of the optical/near-infrared observations, some notes on the cluster properties, 
and original references. We note that some of the target clusters, in particular the highest 
redshift ones, are the targets of wide-field studies by other groups, including the PISCES 
Survey (Kodama et al. 2005). We have included them in our survey because the number 
of well-studied clusters at these redshifts is still limited and because there exists a wealth 
of publically available ground (Subaru and CFHT) and space-based (HST, Chandra, XMM, 
and/or Spitzer) observations on these systems. Thus, we can increase our sample size at 
very little observational cost to us. 

The cluster sample also includes CI 1604+4304 at z = 0.90 and CI 1604+4321 at 
z = 0.92, the original clusters detected in the CI 1604 supercluster (see Gal et al. 2008 
and references therein). As the first large scale structure detected in the ORELSE survey, 
the CI 1604 supercluster currently has the most extensive multi-wavelength photometric and 
multi-object spectroscopic data. As discussed below, we use the properties of the 419 con- 
firmed members in the CI 1604 supercluster as a guide when choosing the appropriate galaxy 
magnitude and color selections in our other cluster fields. 

2.2. The Observations 

To detect large scale structure, we need to probe scales greater than 10 Mpc, the 
mean projected separation of cluster pairs in local superclusters (Bahcall et al. 1988). Thus, 
the first phase of our survey is wide field (25') optical/near-infrared (r'i'z'K s ) imaging of the 
20 cluster fields. All of the near-infrared observations are original to the ORELSE survey, 
while the optical imaging comes from both new and archived observations (see Table [T]). 
For our original observations, we are using the Large Format Camera (LFC) and the Wide- 
Field Infrared Camera (WIRC) on the Palomar 5-m, Suprime-Cam on the Subaru 8-m, the 
Wide-Field Infrared Camera (WFCAM) on the UKIRT 3.8-m, and the Florida Multi-object 
Imaging Near-IR Grism Observational Spectrometer (FLAMINGOS) on the KPNO 4-m to 
cover an area of at least 0.2 deg 2 around each target cluster. The optical imaging is initially 
used to characterize the large scale environment around the target clusters and choose prob- 
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able cluster members for follow-up multi-object spectroscopy. The near-infrared K s imaging 
will be used to measure stellar masses and provide improved photometric redshifts. The 
photometric observations reach typical depths of {r' i' V K s } = {24.7 24.2 23.4 20.5} for a 
3a detection. 

Based on the optical (r'i'z') imaging for each cluster field, we use broad color-color 
cuts, characteristic of red, elliptical-like galaxies at the cluster redshift (see §2.3[) . to produce 
adaptively smoothed galaxy density maps. These maps are used primarily to provide a 
visual locator for companion clusters and intermediate-density large scale structures in the 
imaging field (see §2.41) . Filaments and cluster infall regions can cover significant portions of 
the observed area, but at relatively low contrast, making it difficult to detect them and define 
their boundaries. Identification of such overdense regions is used to guide the placement of 
follow-up observations, especially slitmasks for multi-object spectroscopy. 

Because photometric techniques for studying galaxy and cluster evolution (such as pho- 
tometric redshifts, statistical overdensities, and richness estimates) are highly uncertain, 
the crucial component of this survey, and what distinguishes it from competing studies, is 
the unprecedented amount of spectroscopic coverage. We utilize the wide-field, multi-object 
spectroscopic capabilities of the DEep Multi-Object Imaging Spectrograph (DEIMOS; Faber 
et al. 2003) on the Keck 10-m. Because of its large FOV (16.7' x 5.0'), capability of position- 
ing up to 120+ galaxies per slitmask, and high efficiency, DEIMOS is the ideal instrument 
for a spectroscopic survey of high-redshift large scale structures. We perform faint object 
spectroscopy (for objects down to i' « 24) with total exposure times of 2.5-3 hrs. Using 
the density map as a guide (see above) and a series of appropriate color cuts to preferen- 
tially select galaxies at the cluster redshift (see §2.3.21) . we are obtaining high resolution 
(~ 60 km s _1 ) spectra for 100-200+ confirmed members per system. Our sample sizes are 
larger than most of those currently available (e.g., Tanaka et al. 2006, 2007b; Koyama et al. 
2007; Swinbank et al. 2007; Fassbender et al. 2008) and allow us to measure properties on 
a galaxy-to-galaxy basis. Based on the spectroscopy completed so far (see §3 and Gal et al. 
2008), we find that 25-45% of the galaxies targeted through our system of color selections 
are at or near the cluster redshift. 

Based on our imaging and spectroscopic observations, promising supercluster candidates 
and other complex systems are currently being targeted for multi-wavelength studies with 
HST, Chandra, Spitzer, and/or the VLA. We note that our large spectroscopic databases 
are important to capitalize on the multi-wavelength data as we can begin to identify and 
study rare galaxy populations (like X-ray and radio sources) on an individual rather than a 
statistical basis (see Kocevski et al. 2008a,b). 
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2.3. The Color Selection 

To characterize the large scale structure around our target clusters, we need to identify 
galaxies (which will be used to create a density map and which will be targeted for spec- 
troscopy; see below) within the wide-field imaging that are at or near the cluster redshift. 
To do this, we target red, elliptical-like galaxies which provide the best contrast against the 
large background of blue field galaxies. Until we have sufficient spectroscopy in each field 
to refine the color selection, we need to choose initial color cuts in the r'i'z' data which 
are appropriate for the cluster redshift. We have exami ned the possibility o f appl ying color 



cuts based on the stellar population synthesis models of iBruzual fc Charlotl (120031 . hereafter 
BC03). To see the variation in colors between different star formation histories, we generate 
synthetic galaxy spectra using three solar- metallicity models : (1) an instantaneous burst 
with a formation redshift of Zf = 3; (2) a r = 0.6 Gyr exponential decline at Zf = 3; and 
(3) a r = 0.6 Gyr exponential decline at Zf = 5. In Gal et al. (2008), we also exp lored the 



prescription used by the Red-sequence Cluster Survey (RCS. lGladders h Yedl2005l ). a model 
parameterized as a 0.1 Gyr burst ending at z — 2.5, followed by a r = 0.1 Gyr exponential 
decline. We do not examine the RCS model here as it is nearly identical (to within 0.04 mag) 
to the simpler model (#1) for the redshift range in which we are interested (0.6 < z < 1.3). 

All three models reproduce the colors of local elliptical galaxies (e.g., Blanton et al. 
2003; Chang et al. 2006) and have formation epochs implied by the observed color evolution 
in red-sequence galaxies (e.g., Stanford et al. 1998; van Dokkum & Franx 2001; Blakeslee 
et al. 2006; Mei et al. 2006b; Homeier et al. 2006). The evolution with redshift in the 
(r' — i') vs. (i' — z') color-color diagram (CCD) is shown in Figure [1] for each of the three 
models. The first thing to note is that all three models are identical at z < 0.6 and nearly 
identical at redshifts of 0.6 < z < 1, suggesting that the color selection at these redshifts is 
insensitive to the model as long as the formation epoch is sufficiently high. We do, however, 
observe significant variations at z > 1, clearly due to differences in formation redshift and 
starburst length. We note that the abrupt changes in color with redshift, most notably in the 
instantaneous burst model (#1), are due to strong Fell and Mgll absorption lines shifting 
in and out of the optical bandpasses. 

Overlaid in Figure [1] are the average colors of spectroscopically confirmed red-sequence 
galaxies in three ORELSE structures, the CI 1604 supercluster at z ~ 0.9 (Gal et al. 2008), 
the CI 1324 supercluster at z « 0.73 (Gal et al. in prep), and RX J1821.6+6827 at z = 0.82 
(see Figure [3]). We do not include CI 0023 here as it is dominated by blue galaxies and 
has a wider red sequence than the other systems (see §4. 2D . This figure demonstrates that 
a single model does not fit all three systems. In particular, the colors from the CI 1604 
supercluster and RX J1821. 6+6827 are redder in r' — i' and bluer in i' — z' than all of the 
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models. Such discrepancies are unsurprising as the evolutionary history of these systems 
are undoubtedly different and more complex. In addition, a galaxy's star formation history 
will vary with environment. As a result, we adopt broad (Am = 0.4) color cuts that were 
used successfully in the CI 1604 supercluster. We use the models as a guide and offset our 
color cuts relative to those used in CI 1604 (see §2.2 of Gal et al. 2008). For cluster redshifts 
between z = 0.6 — 1.0, we use a standard r' — i' cut of 1—1.4 but vary the i' — z' cut depending 
on the cluster redshift. As done for the CI 1604 supercluster, we adopt an i' — z 1 cut of 0.6-1.0 
for a cluster at z = 0.9 — 1.0. For those clusters at lower redshifts, we adopt an i' — z' cut 
which gets progressively bluer by 0.1 mag as the redshift interval decreases by Az = 0.1. 
For cluster redshifts of z > 1, we use a fixed cut of i! — z' > 0.9. We use a single color cut 
at these redshifts because, while the r' — i' color at z > 1 is highly dependent on the model, 
the i' — z' colors are consistently above ~ 0.9, except for models with recent or continuing 
star formation. We note that ACS observations of three X-ray-selected clusters at z > 1 all 
reveal well-defined red sequences with characteristic colors of zf775W — ^fssolp ~ 1 (Mei et 
al. 2006a,b; Demarco et al. 2007). These results are completely consistent with our adopted 
%' — z' color cut. Table [2] lists the exact color cuts used in this survey. 

We note that our red-galaxy selection plus density-mapping technique is fundamen- 
tally differ ent from cluster det ection methods such as the red-sequence (RCS) method and 



MaxBCG (IKoester et al.l 120071 ). We use a simple, relatively wide (Am > 0.4 mag) color cut 
to examine the projected density of objects with colors broadly consistent with those of red 
galaxies at the redshift of interest. RCS and MaxBCG both rely on the presence of a tight 
red sequence (with an assumed intrinsic scatter of 0.075 and 0.060, respectively), which, 
especially at higher redshift, may impose a priori requirements on the types of systems de- 
tected. In the ORELSE survey, we target fields with known structures and search for lower 
significance density enhancements, such as groups and filaments, whose constituent galaxies, 
while more evolved (redder) on average than field galaxies, do not necessarily form a tight 
red sequence (as is the case for CI 0023; see §4.2p . While we have made broad assumptions 
on the formation history of the cluster galaxies when choosing our color cuts (particularly at 
z > 1) in order to eliminate fore/background galaxies, the redshifts of our targets are already 
known, so we do not require a photometric redshift estimator such as the red-sequence fitting 
(which relies on one particular evolutionary model) used by RCS. 



2.3.1. Magnitude Limits for the Density Map 



In addition to the color cuts used to select red galaxies at the cluster redshift, we also 
need to define the range of il magnitudes used to produce the density map (see below). In 
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the case of the CI 1604 supercluster at z — 0.9, we chose %' = 20.5 — 23.5 based on our 
existing spectroscopy. For consistency, we shift the %' magnitude limits to the redshift of the 
other target clusters by using /c-corrections derived from the instantaneous burst at zj = 3 
model (#1) shown in Figured) The magnitude limits used to create the red-galaxy density 
maps for the two target clusters presented in this paper are given in §3.31 



2.3.2. Priorities for Spectroscopic Observations 

When determining which galaxies to target for multi-object spectroscopy, we prioritize 
galaxies based on color. As a first priority, we select objects which meet our color cut(s) for 
red galaxies described above. Within this priority (and the other lower priorities), we weight 
objects by their i'-band flux with preference given to brighter objects. We include objects 
with i! between the minimum magnitude used to create the density map (see §2.3. ip and 
i' = 24, the limit where we can still obtain reasonably high signal-to-noise (S/N ~ 1 — 5) 
spectra in our 2.5-3 hr exposures. With increasingly lower priority, we select galaxies which 
are progressively bluer in both r' — i' and i' — z' to target galaxies within the green valley 
and into the blue cloud. Although our contamination rates are higher for bluer galaxies, we 
observe a significant population of blue, star-forming members in the CI 1604 supercluster 
(see Figure 1 of Gal et al. 2008) and CI 0023 (see Figure |2j), so it is essential to target these 
galaxies as well. We note that we always maximize the number of slits per slitmask (given 
our minimum slit length of 5") by including all other galaxies, within our magnitude limits, 
at the lowest priority. The placement of these slitmasks (with the prioritized galaxy selection 
described above) on the sky is guided by the structure (i.e., overdense regions and filaments) 
identified in the red-galaxy density map (see §2.41) . 



2.4. Producing the Density Map 



Following iGal fc Lubinl (120041 ). we produce a red-galaxy density map for each cluster 
field by adaptively smoothing a subset of galaxies which satisfy the galaxy color selection 
and magnitude limits determined in §2.31 An adaptive kernel smoothing is applied using an 
initial window of 0.75 h^ Mpc and 10 arcsecond pixels. The chosen kernel size prevents 
small groups from being blended into single detections in the density map and enhances the 
contrast of small groups against the background, making detection of such low-mass systems 
easier. We stress that the density maps of color-color selected galaxies are primarily a tool 
for the visual identification of structures such as clusters, groups, walls, and filaments. 
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2-4- 1- The Cluster Detection Threshold 



Because the ORELSE survey aims to detect comparatively low mass systems, our clus- 
ter/group selection algorithm is potentially susceptible to a high false detection rate. It is, 
therefore, useful to select overdensities on the basis of some objective criteria. The applica- 
tion of both r' — i' and i' — z' color cuts to select galaxies for the density maps results in a 
non-trivial relationship between projected galaxy density and the cluster mass. A complete 
treatment of selection functions would require large cosmological simulations with model 
galaxy colors in our photometric system and would still depend on the assumed cosmology 
and recipes for modeling galaxy formation. Therefore, we rely instead on two observational 
datasets: t he extensive photometry and spectro scopy in the field of the CI 1604 supercluster 
at z = 0.9 ( Gal fc Lubinl 12004 iGal et al.ll2008l ) for the cluster component and the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) Deep Fields^ for the "field" component. 

The CI 1604 system consists of numerous groups and clusters, with velocity dispersions 
as low as ~ 300 km s _1 , corresponding to a mass of ~ 10 14 hy M Q . The spectroscopic 
catalog of over 400 confirmed supercluster members allows us to reliably identify such low 
mass systems and examine their prope r ties in the map of galaxy density. We use the structure 
identified as Cluster C in IGal et al.l (120081 ) to set our detection thresholds because it is 
the lowest-mass group to be spectroscopically confirmed in the CI 1604 supercluster. This 
poor cluster is evident as a single peak in the density map and has a sufficient number of 
spectroscopic me mbers (18 within 1 h jn Mpc) to estimate its velocity dispersion o~c = 313 
km s _1 . Following IGal &: Lubinl (12004 ). we use SExtractor to detect structures in the density 
maps. The SExtractor parameters must be chosen with care, to recover systems which we 
know to be real (in this case, Cluster C and all of the more massive systems in CI 1604), 
while minimizing the number of spurious detections. 

To measure the false detection rate, we require a galaxy distribution that retains the 
general correlation properties of galaxies meeting our color-color selection criteria. Optimally, 
we would measure the correlation function of field galaxies from the same data used for 
cluster selection. However, our imaging survey is entirely targeted at known high-density 
regions of the universe. Thus, we use the CFHTLS to derive the statistical properties of the 
red galaxy distribution. The 2008A data release includes four fields, of ~ 2 square degrees 



1 Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope (CFHT) which is operated by the National Re- 
search Council (NRC) of Canada, the Institut National des Science de l'Univers of the Centre National de la 
Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data 
products produced at TERAPIX and the Canadian Astronomy Data Centre as part of the Canada-France- 
Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. 
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each, at widely separately locations on the sky. Having four distinct fields allows us to 
estimate not only the typical false detection rate but also the range due to cosmic variance. 
Matched photometric catalogs for all four fields are provided, including r^j/z' p hotom etry 



(the same as our d ata) in the AB system. In brief, following iPostman et al.l (119961 ) and 



Gal fc Lubinl (120041 ). we use the CFHTLS data to compute Raleigh-Levy (RL) parameters 
for galaxies meeting our color and magnitude cuts which are then used to generate simulated 
galaxy distributions. These distributions are used to produce density maps, on which we 
run SExtractor. The SExtractor parameters DETECT_THRESH and MIN_AREA are varied 
until a set of values is found that successfully detects the confirmed groups/clusters in CI 
1604, while yielding low contamination. 

In each CFHTLS Deep Field, we select those objects classified as galaxies, meeting the 
same quality, i' magnitude, and r' — i', i' — z' cuts applied to the Palomar LFC data for 
the CI 1604 imaging. This gives the number of galaxies Nrlj that should be used in the 
RL simulation for the /-th CFHTLS Deep Field. This number varies for every ORELSE 
target as the magnitude and color limits change. A square region of exactly 1° on a side 
is selected from the center of each CFHTLS Deep Field to simplify computation of the RL 
parameters and avoid underexposed field edges. We then generate a set of simulated RL 
distributions over a grid of values for the two RL parameters 9q and d. The best match to 
the observed CFHTLS data yields values of #0 = 100", d = 0.6 using the color and magnitude 



cuts from CI 1604. The se parameters are comparable to those found by iGal fc Lubinl (120041 ) 



and [Lopes et ajj (120061 ) 



The next step is to create 100 simulated galaxy distributions for each Deep Field using 
these RL parameters, each simulation having Nrlj galaxies. This yields 400 total density 
maps. We then run the same adaptive kernel density mapper applied to the real data on 
these simulated distributions and detect overdensities using SExtractor. We vary the DE- 
TECT_THRESH and DETECT_MINAREA parameters, using values of the galaxy density 
at different radii r test from the center of Cluster C in CI 1604, setting DETECT_MINAREA 
to 7rrf est . We find that an area of 47 pix 2 , corresponding to a radius of 0.3 Mpc, allows 
us to detect all of the spectroscopically confirmed cluster candidates in the real data, while 
minimizing the false candidates detected in the RL simulations. The corresponding galaxy 
density threshold is 7300 galaxies per square degree. The median number of false detections 
expected in the area imaged for CI 1604 is Nf a i se = 0.6, compared to the ten candidates 
detected (see Figures 2-3 of Gal et al. 2008). The range of false detections based on the four 
CFHTLS Deep Fields is between 0.3 and 0.9 for the CI 1604 area. 



As demonstrated in ILopes et al.l (120061 ) . the RL simulations provide a good estimate 



of the false detection rate. A much simpler approach is to generate purely random galaxy 
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distributions, but these do not encode any of the galaxy correlations present in the real 
universe. Nevertheless, they serve as a useful confirmation of the RL results. We performed 
the same analysis described above, using random positions instead of an RL distribution. 
As expected, we detect many fewer false candidates, with Nf a i se = 0.4 for the CI 1604 area. 

2.4-2. False Detection Rates 

Having set the SExtractor detection parameters based on CI 1604, we run cluster detec- 
tion on the density maps for other cluster fields with the same parameters. However, because 
the color and magnitude limits vary from field to field based on the cluster redshifts and 
pre-existing spectroscopy, the RL parameters and especially the number of galaxies used in 
the RL simulation for each field will vary, resulting in different contamination rates. First, 
we followed the same procedure described above to determine the optimal RL parameters 
for each ORELSE targets. For all such targets, we find that the range of well-matched pa- 
rameters is quite broad, so we fix the RL parameters to 9 = 100", d = 0.6 as found for 
CI 1604. Therefore, the only quantity that we change in our RL simulation from target to 
target is the total number of galaxies. 

To determine the number of galaxies N RL to use for each ORELSE target, we simply 
apply the %' magnitude limits and r' — i' and i' — z' color cuts used for the target to each of 
the four CFHTLS Deep Fields. 

For every target field t, we proceed as follows: 

1. We generate 100 RL distributions over each square CFHTLS Deep Field with N RLjt 
galaxies. This number changes between the Deep Fields due to cosmic variance. 

2. Density maps are created from each of the distributions, with the same 0.75 Mpc 
smoothing used for the corresponding target. 

3. SExtractor is run with DETECTJvlINAREA corresponding to a radius of 0.3 Mpc 
at the redshift z t of the target, and DETECT_THRESH fixed at the galaxy density 
from Cluster C in CI 1604. 

4. The number of expected false detections Nf a i sejt is set to the mean value from the 400 
simulations, scaled to the area imaged for target t. 

The median number of false detections (and their range from the four CHFTLS Deep 
Fields) expected in the fields of the two target clusters discussed in this paper, CI 0023 at 
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z = 0.84 and RX J1821.6 at z = 0.82, is 0.60 (0.3-1.0) and 0.55 (0.3-0.8), respectively. 
Cosmic variance on the scale of the LFC field could still potentially cause significant changes 
in the number of false detections. We, therefore, repeated the above process, now using 100 
randomly placed LFC-sized subfields within each of the 400 RL density maps for each target. 
We find that the median and rms number of false detections is 0.5 ± 0.8 and 0.3 ± 0.6 for CI 
0023 and RX J1821, respectively. However, the median richnesses of these false detections 
is typically much lower than for the overdensities found in our target fields (see §1]). 



3. The Targets and Observations 

3.1. CI 0023+0423 

CI 0023 at z = 0.84 was originally detected as an overdensity in the distant cluster survey 
of Gunn, Hoessel & Oke (1984) which covered 71.5 deg 2 using photographic plates taken with 
the Palomar 1.2-m and 5-m and the KPNO 4-m telescopes. This system was later included in 
the Oke, Postman & Lubin (1998) survey using the Low Resolution Imaging Spectrograph 
(LRIS; Oke et al. 1995) on the Keck 10-m. This study spectroscopically confirmed 24 
members and revealed that this system is actually comprised of two small groups, with 
velocity dispersions of 158 and 415 km s _1 and separated by 2922 km s _1 in radial velocity 
(Postman, Lubin & Oke 1998). N-body simulations implied that the two groups were likely 
bound, were curre ntly merging, and w ould form a massive cluster (ai os ~ 730 km s _1 ) in 
the next ~ 1 Gyr flLubin et al.lll998bl ). Based on the LRIS photometry and spectroscopy, 
Postman, Lubin & Oke (1998) found that the population is quite active, with 57% of the 
confirmed members showing [Oil] emission with equivalent widths greater than 15 A and 
the majority with broad-band spectral energy distributions ("color ages") which suggest 
star formation within the past 3 Gyr. Twelve confirmed members were covered by a single 
WFPC2 (F702W) image centered on this system. 75% are late-type galaxies, a morphological 
fraction more consis tent with galaxy groups or the field than rich clusters at these redshifts 
flLubin et al.1ll998ah . 



Follow-up multi-wavelength studies of this system also indicate substantial star forma- 
tion and AGN activity. SCUBA observations of four high-redshift (0.84 < z < 1.27) clusters, 
including CI 0023, reveal 10 securely-identified submm sources detected toward the cluster 
fields, a combined source excess at 850 /im of ~ 3 — 4 times blank field surveys (Best 2002). 
Similarly, an H Q imaging survey of CI 0023 revealed a high integrated star formation rate 
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observed in H Q cluster surveys at z ~ C 
not find a population of early-type galaxies in CI 0023 with on-going star formation (i.e. 
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detectable H a emission), suggesting that the [Oil] emission that we observe in otherwise red, 
passive galaxies may, in fact, be due to AGN activity (see §4.21) . 



3.2. RX J1821.6+6827 

RX J1821 at z = 0.82 was the highest-redshift cluster discovered in the ROSAT North 
Ecliptic Pole (NEP) survey (Gioia et al. 2003; Henry et al. 2006). XMM observations revealed 
slightly elongated X-ray emission, with a bolometric luminosity of 1.17lo']J x 10 45 /if 2 erg 
s _1 and a temperature of 4.7lo!7 keV (Gioia et al. 2004). Spectroscopic redshifts for twenty 
cluster galaxies using the Keck I and the Canada-France-Hawai'i (CFH) telescopes imply a 
cluster velocity dispersion of 775^}^ km s _1 , typical of a relatively rich cluster (Gioia et al. 
2004). 

The cores of all 18 clusters in the NEP survey were imaged by the VLA at 1.4 GHz 
(Branchesi et al. 2006). A total of 32 radio sources, with peak brightnesses > 0.17 mJy/beam, 
were detected, suggesting a 3a overdensity of sources relative to the field within 0.4 Mpc. 
One of these sources was a confirmed member of RX J1821. These results indicate radio 
AGN activity in the cluster cores, consistent with studies of other clusters (Owen et al. 1999; 
Best 2003; Best et al. 2002; Miller & Owen 2003). 



3.3. Photometric Data 



CI 0023 and RX J1821 were imaged using the Large Format Camera (LFC; ISimcoe et al. 



20001 ) on the Palomar 5-m telescope. The LFC is a mosaic of six 2048 x 4096 CCDs in a 
cross-shaped layout, mounted at prime focus. The imaging area corresponds roughly to the 
unvignetted region of a circle 24' in diameter. Data were taken on UT 27 April and 15 
August 2004, with the camera installed in the standard North-South orientation and all six 
CCDs in use. The imaging data were taken in unbinned mode, with a pixel scale of 0'/182 
pix -1 . Average seeing was l'.'O, and the nights were not photometric. Individual exposures 
of 450 sec were taken, moving the telescope by up to 45" in each direction to alleviate the 
gaps due to the spaces between the CCDs. Each target was imaged in the Sloan Digital Sky 
Survey (SDSS) r'i'z' filters, with total exposure times in these filters of {4500 3600 2250} 
sec for CI 0023 and {6300 4050 3600} sec for RX J1821. 



Data reduction was performed using the IRAF data reduction suite ( ITodylll986l ). includ- 
ing the external package mscred. We followed the general mosaic imaging guidelines used 
by the NOAO Deep Wide Field Survey (IJannuzi et al.ll2004l ). with modifications appropriate 
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to our data. The details of the data reduction are described in iGal et al.1 (120051 ). Object 
detection in the target frames was performed by running SExtractor in dual-image mode, 
using ultra-deep images (combining all three bands) for detection, while measurements were 
performed on the single-band images. Elliptical photometry apertures are determined from 
the deep image, and these same apertures are used for the individual filter images. This 
procedure results in catalogs where magnitudes are measured using identical apertures in all 
three filters, which improves the measurement of g alaxy colors by using the same physical 



size for each galaxy in all filters (ILubin et al.1 120001 ) . Due to the lack of calibration, we as 



sumed photometric zero points, color terms, and atmospheric extinction terms from other 
photometric nights, resulting in semi-calibrated catalogs. 

To provide photometric calibration, RX J 1821 was observe d with the Orthogonal Par- 



allel Transfer Imaging Camera (OPTIC; iHowell fc Tonryl 120031 ) installed on the University 



of Hawai'i 2.2-m telescope on UT 22 July 2007. This instrument contains two 2K x 4K 
CCDs mounted side-by-side, with 0'/14 pixels. While this instrument is capable of orthog- 
onal transfer charge shifting to improve image quality, we used it in conventional mode to 
minimize data reduction complexity. Three exposures of 300 sec each were obtained in each 
of the SDSS r', i', and z' filters in photometric conditions with 0'.'8 seeing. Numerous SDSS 
standard stars were also observed to derive photometric calibration. These observations 
were processed and sources detected the same way as for our LFC data. The calibrated, 
photometric source catalogs from OPTIC were matched to the semi-calibrated catalogs from 
LFC, and transformations of the form 

m ca ii b = m LFC + A x (color) + B (1) 

were derived for each filter, using the r' — i' color for the r' and i' filters, and i' — z' color for 
z' . These transformations were then applied to the full source catalogs from the LFC data. 
Unfortunately, we have not yet been able to obtain photometric calibration for the CI 0023 
field. Comparison of the semi-calibrated LFC data of RX J1821 with the final calibration 
using OPTIC suggests that zero points in each filter may be incorrect by 0.05-0.1 mag, with 
errors of ~ 0.1 on the color term. 

The resulting color-magnitude diagrams (CMDs) and color-color diagrams (CCDs) for 
the CI 0023 and RX J1821 fields are shown in Figures [2] and [3j Based on the photometry 
described above, red-galaxy density maps were made using the appropriate color-color cuts 
(see §231 and Table EJ and i'-band magnitude ranges of 20.28-23.28 and 20.16-23.16 for CI 
0023 and RX J1821, respectively (see §2.3.ip . Figures @] and show the resulting density 
maps. Using the RL simulation of §2.41 we expect 0.60 and 0.55 false detections in the CI 
0023 and RX J1821 fields, respectively, significantly fewer the actual number that we detect 
(see @. 
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3.4. Spectroscopic Data 

Both CI 0023 a nd RX J 1821 wer e observed with the Deep Imaging Multi-object Spec- 



trograph (DEIMOSi lFaber et al.ll2003l ) on the Keck 10-m telescopes. We used the 1200 1/mm 
grating, blazed at 7500 A, with the OG550 filter and 1" slits. This setup results in a pixel 
scale of 0.33 A pix -1 , a resolution of ~1.7 A (68 km s _1 ), and typical spectral coverage from 
6250 A to 8900 A. Galaxies as faint as i' ~ 24 were observed, with highest priority given 
to red galaxies based on the color cuts used to produce the initial density maps from the 
semi-calibrated photometry. Approximately 60% of the targeted objects met the red galaxy 
color cut. Lower priorities were given to galaxies increasingly blueward (in r' — i' and/or 
i' — z') of the highest priority box. CI 0023 was observed with three masks and RX J1821 
with two masks on UT 02-03 September 2005 in / .'7 seeing, with 5 exposures of 1800s each 
on each mask. CI 0023 was observed with two additional masks on UT 12 September 2007 
with 5 exposures of 1800s each in 0'.'8 seeing and thin cirrus. 

Data were reduced usi ng the DEEP2 version of the spec2d and specld data reduction 



pipelines (IDavis et all 120031 ) . Redshifts were obtained using DEEP2's redshift measurement 
pipeline zspec (Cooper et al. in prep), with typical redshift errors of 25 km s -1 . Each redshift 
was assigned a quality flag between 1 and 4, where 1 indicates a secure redshift could not 
be determined due to poor signal, lack of features, or reduction artifacts; 2 is a redshift 
obtained from either a single feature or two marginally detected features; 3 is a redshift 
derived from at least one secure and one marginal feature; and 4 is assigned to spectra 
with redshifts obtained from several high signal-to-noise features. A quality of -1 is used for 
sources securely identified as stars. Qualities of 3 or 4 indicate secure extragalactic redshifts 
and only those are used in the following analyses (see also Gal et al. 2008). 

In the CI 0023 field, we obtained 423 extragalactic redshifts with Q > 3, of which 54 
were serendipitous detections. An additional 73 galaxy redshifts were incorporated from the 
spectroscopic survey of Oke, Postman & Lubin (1998). The combined catalog contains Q > 3 
redshifts for 169 structure members with 0.82 < z < 0.87. In the RX J1821 field, a total 
of 129 extragalactic redshifts with Q > 3 were obtained in the RX J1821 field, 26 of which 



were serendipitous detections. Thirteen additional redshifts were taken from iGioia et al. 



(120041 ) . The combined catalog contains Q > 3 redshifts for 85 structure members with 



0.80 < z < 0.84. 
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4. The Results 

4.1. The Global Properties 

In Figures H] and we show the red-galaxy density maps for the CI 0023 and RX J1821 
fields. The density peaks which meet our detection threshold (see §2.4. 1[) are circled. We find 
three peaks in the CI 0023 field and one in the RX J1821 field, all significant based on our 
estimated false detection rate. The observed structures in these maps are consistent with 
that expected from the redshift distributions shown in Figures El and [7J While Postman, 
Lubin & Oke (1998) found only two peaks in the redshift distribution, we see four clear 
peaks with our more extensive spectroscopy. Combined with the density map, these results 
suggest an even more complex system of groups. In contrast, the redshift histogram of RX 
J1821 shows a single peak at z ~ 0.82 (albeit with some clear velocity substructure; see 
below) which is consistent with the single, massive overdensity observed in the density map. 

In addition, the distribution of confirmed members in each system follows exceedingly 
well both the peaks and lower-density regions in the two density maps. In the CI 0023 field, 
although many of the confirmed members are spread out across the spectroscopic area, 54% 
coincide directly (within 1 Mpc) with the three significant peaks (labeled A, B, and C). 
This correlation is not due to spectroscopic coverage alone as the five DEIMOS masks cover 
a large fraction of the imaging area, with most of that area covered by at most two masks 
(see left panel of Figure HJ). We also see that the density map does an excellent job at picking 
up associated lower-density structures, such as the less-significant peak directly to the South 
of component B. Splitting the confirmed members by redshift according to the four redshift 
peaks in Figure [6l we find velocity segregation between the three components. Component A 
is largely comprised of galaxies belonging to the redshift peak at z ~ 0.839, while component 
C is largely comprised of galaxies belonging to the redshift peak at z ~ 0.845. Component 
B, however, is actually a superposition of two groups along the line-of-sight, one which 
corresponds to the lowest redshift peak at z pa 0.828 and the other which corresponds to the 
redshift peak at z pa 0.845, same as component C. The galaxies within the highest redshift 
peak at z ~ 0.864 are not centrally concentrated but rather extend across the entire region, 
suggesting a sheet of galaxies in the near background. 

In Figure [8], we show the redshift histograms, within radii of 0.5 and 1.0 Mpc, 
centered on the three overdense components. Following the procedure described in Gal et al. 
(2008), we use these distributions to measu re the biweight velocity dispersion computed by 



ROSTAT (IBeers. Flynn. Sz Gebhardtlll990l ). Components A and C represent single systems 
with measured velocity dispersions between 428-497 km s _1 , typical of massive groups or 
poor clusters (e.g., Zabludoff & Mulchaey 1998). Component B is comprised of two smaller 
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galaxy groups with measured velocity dispersion of 206-293 km s _1 (see Table [3]). As such, 
CI 0023 is a four-way group-group merger. Based on virial mass estimates of the individual 
groups, the final mass of the merged cluster will be ~ 5 x 10 14 M Q , typical of a rich 
cluster (e.g., Girardi et al. 1998). 

While CI 0023 is dynamically active, RX J1821 is considerably more relaxed with a single 
peak in both the density map and redshift histogram (see Figure [5] and [7]) ; however, there 
are indications of continuing formation in this cluster as well. In particular, the overdense 
region in the density map corresponding to the cluster is clearly extended in the North-South 
direction, consistent with the X-ray contours from Gioia et al. (2004). As in CI 0023, the 
distribution of confirmed members follows extremely well the overdense regions in the density 
map, tracing the elongated structure of the cluster and even picking up the lower-density 
peaks (small groups) to the South and South-west which are kinematically associated with 
the cluster (right panel of Figure [5]) . 

To determine if RX J1821 has any velocity segregation in the plane of the sky, we 
have binned the velocities by sky position going from North to South along the elongated 
structure. We use four bins, each with an equal number of galaxies. Comparing each resulting 
velocity histogram to the full velocity histogram using the KS test, only the Southern-most 
distribution is marginally different, with a 6% probability of being drawn from the same 
distribution, possible due to an infalling group (see Figure [5]). Otherwise, no statistically 
significant velocity segregation was detected. 

Velocity histograms, centered on the density peak and measured within a projected 
radius of 0.5, 1.0, and 1.5 Mpc, are shown in Figure [9j We measure biweight velocity 
dispersions between 926-993 km s -1 (see Table |3]). We note that three ROSTAT dispersion 
estimators (sigma, biweight, and gapper) give consistent measurements within a radial bin 
and from bin to bin; however, as is clear from Figure [HJ the velocity distributions are clearly 
non-gaussian, implying measurable velocity substructure. RX J1821 does not, however, have 
clear multiple components like CI 0023 or a significant amount of nearby structure like the 
CI 1604 supercluster (Gal et al. 2008). 

The variation in system properties might reflect the differences in selection technique. 
Based on the handful of ORELSE targets observed so far, we find that optical selection 
appears to preferentially select systems which are dynamically active (CI 0023), are embedded 
in complex, large scale structures (the CI 1324 and 1604 superclusters) , contain lower-mass 
constituents (all of the above), and have less (or no) diffuse gas emission (all of the above; 
see Lubin, Mulchaey, & Postman 2004; Kocevski et al. 2008a; Kocevski et al. in prep). As 
often suggested, X-ray selection appears to detect more dynamically-evolved, concentrated, 
and isolated clusters (RX J1821 and another NEP cluster, RX J1757.3+6631 at z = 0.691) 
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at these redshifts. 



4.2. The Galaxy Populations 

Although CI 0023 and RX J1821 are at almost identical redshifts, they are distinctly 
different in their evolutionary state. As a four-way group merger, CI 0023 is at the very early 
stage of hierarchical cluster formation. In comparison, the cluster RX J1821 has progressed 
further along its dynamical evolutionary track as evidenced by its more concentrated velocity 
and spatial structure. These global differences translate to substantial differences in their 
galaxy populations as well. Based on the CMDs in Figures [2H3] and the color histograms 
in Figure [T0| CI 0023 has a significantly larger population of blue galaxies, with 51% being 
bluer than the red-galaxy color-color cut (down to i' = 24.5), compared to only 24% in RX 
J1821. RX J1821 also has a tight, well-defined red-sequence which extends down to at least 
i' ~ 23.5, unlike the wider one of CI 0023 (although its larger width may be due, in part, 
to photometric errors; see §3.31) . The CMD of RX J1821 indicates a population of both 
old, massive (bright) red galaxies with high formation epochs (zf ~ 2 — 3), consistent with 
previous studies of massive clusters (e.g., Mei et al. 2006b and references therein), as well as 
quenched galaxies as evidenced by the presence of faint red-sequence galaxies (e.g., Koyama 
et al. 2007 and references herein) and stronger H<5 absorption (see below). 

The diversity in activity is also evident in the galaxy spectra. In CI 0023, 80% of the 
confirmed members have detectable [Oil] emission (equivalent widths greater than 2 A), 
including detectable line emission in otherwise red, passive galaxies (with Ca H&K but no 
Balmer absorption). We note that some of the [Oil] emission, in particular in these passive 
galaxies, may actually be due to AGN and/or extended emission-line regions, i.e., Seyferts 
or LINERs (e.g., Yan et al. 2006). This high [Oil] fraction is consistent with the results 
of Poggianti et al. (2006) who studied a sample of groups and clusters at z — 0.4 — 0.8 
from the ESO Distant Cluster Survey (EDisCS). They find an anti-correlation between [Oil] 
fraction (/roiii) and velocity dispersion (er), with the majority of groups (a < 500 km s _1 ) 
having /[on] > 0.7 (see Figure 4 of Poggianti et al. 2006). Thus, the high [Oil] fraction in CI 
0023 could just result from the fact that the constituents of the system are all group-sized. 
In contrast, the population of the massive cluster RX J1821 is more quiescent, with only 
36% of the galaxies showing [Oil] emission. The [Oil] fraction is completely consistent with 
that found in EDisCS clusters with similar velocity dispersions. Although comparable to 
other moderate-to-high redshift clusters, we note that the measured fraction in RX J1821 is 
higher than the /[on] ~ 0.23 found in nearby (z = 0.04 — 0.08) SDSS clusters with similar 
dispersions (see Figure 6 of Poggianti et al. 2006), consistent with the increase in blue, 
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star- forming cluster galaxies with redshift (e.g., Butcher & Oemler 1984; Ellingson et al. 
2001). 

In Figure [TT], we show a comparison between the composite spectra from the CI 0023 
and RX J1821 galaxy populations. When making the composites, we use only the DEIMOS 
spectra because of their higher quality and spectral resolution. Each spectrum is normalized 
(by their average rest-frame 3700-4300 A flux) and weighted by its z'-band luminosity. A 
small fraction (less than 13%) of galaxies are excluded because of instrumental effects or 
bad columns in their spectra or non-detections in the photometry. The resulting composites 
include 138 and 69 galaxies in CI 0023 and RX J1821, respectively. The comparison reveals 
clear differences between the two galaxy populations. The CI 0023 composite has stronger 
[Oil] and H5 lines, as expected from the overall blue colors of the galaxy population. There 
is also evidence of significant infill of the H<5 line from ongoing star formation which is 
essentially absent in the RX J1821 data. Without this infill (or equivalently ~ 100 Myr after 
star formation is truncated), the apparent contribution from A stars (i.e., the post-starburst 
mode) would be much higher than for RX J1821. Conversely, the composite of RX J1821 
has strong Ca H&K lines and a significant excess of flux at wavelengths greater than 4500 A, 
showing a stronger contribution from K and M stars. Using the standard bandpasses from 
Fisher et al. (1998), we measure the equivalent widths (EWs) of the [Oil] and B.5 line for 
each composite spectrum. We find {OII,H5} = {-8.0,3.9} for CI 0023 and {-4.7, 1.4} for 
RX J1821. Using slightly different bandpasses (from ejg., Postman et al. 1998 and Balogh 
et al. 1999) or the line-fitting functionality of the irafg task splot changes the EWs by less 
than 20% for B5 and less than 10% for [Oil]. 

Because the spectroscopy in CI 0023 is more extensive than that in RX J1821 (i.e., 
5 versus 2 slitmasks, respectively), it is possible that we probe further into the bluer, low 
spectroscopic priority population, thus skewing the comparison of the two systems. To test 
this possibility, we compared only the data from the first two slitmasks for each of the two 
fields. In these subsamples, CI 0023 has 85 confirmed members, 54 of which (64%) meet our 
red galaxy color cut, whereas RX J1821 has 69 confirmed members, with 56 (81%) meeting 
this color cut. These percentages, along with Figures [2H3J confirm the higher red galaxy 
fraction in RX J1821. The excess red galaxies in RX J1821 are faint, red sequence galaxies 
which are much less common in CI 0023, where we instead find more luminous blue galaxies. 
The EWs of the [Oil] and B5 lines in the CI 0023 composite are essentially unchanged from 
the full sample. Even if we consider only the high-priority (red) galaxy populations from the 



2 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for Research in Astronomy under cooperative 
agreement with the National Science Foundation. 
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first two masks, we find {OII,H5} = {-5.8,3.4} for CI 0023 and {-2.9, 1.0} for RX J1821. 
As expected, the overall line strengths are weaker in the red galaxies than in the full sample, 
but the differences between the two clusters persist. Even the red galaxy population in CI 
0023 is more active in star-forming/AGN activity than it is in RX J1821. 

In Figure [T2], we place the results of our composite spectra (from the total galaxy sam- 
ples) on the Dressier et al. (2004) plot of [Oil] vs. H5 equivalent width. The two dashed 
lines, based on the CNOC2 and 2dF field surveys, indicate the average spectral properties 
of a galaxy population which consists of passive galaxies plus a varying fraction of contin- 
uously star-forming galaxies. The low-redshift 2dF data represent the low star formation 
rate (SFR) case, while the moderate-redshift (z ~ 0.4) CNOC2 data represents the high 
SFR case. The small asterisks on the curved lines indicate subsequently higher fractions of 
normal star-forming galaxies, ranging from 20 to 100% of the total cluster population. Local 
clusters from the Dressier & Schectman (1988; hereafter DS) sample fall directly on the lines, 
indicating a cluster population largely dominated by quiescent galaxies plus a small contri- 
bution of continuously star-forming galaxies. Moderate-redshift ({z) = 0.442) clusters from 
the Morphs sample (Dressier et al. 1999), as well as the massive cluster MS 1054 at z — 0.83 
(van Dokkum et al. 2000), all have average H<5 strengths of ~ 2 A and have substantial 
[Oil] emission ranging from -3 to -10 A. Both features indicate significantly more on-going 
star formation than found in the local sample (see Dressier et al. 2004). Unlike the nearby 
clusters, the moderate-redshift clusters lie above the lines, implying an B.5 line which is too 
strong to be produced by normal star formation. These results indicate a larger contribution 
from galaxies that have undergone a recent starburst. 

When placed on this diagram, we see that the average properties of the CI 0023 pop- 
ulation is on the high end of the [Oil] distribution from the Morphs sample and, more 
significantly, substantially offset in B.S, above both the normal star-forming lines as well as 
the Morphs average. These results imply an even larger contribution from recent starbursts 
in the CI 0023 system. In contrast, the RX J1821 population has a weaker [Oil] line, nearer 
to the low-end of the Morphs sample, and falls considerably closer to the normal star-forming 
lines. These results indicate a large quiescent galaxy population, a small population of con- 
tinuously star- forming galaxies (as also evident in the CMD of Figure [3D, and only a modest 
contribution from starburst galaxies. As noted by Dressier et al. (2004), there is strong 
evolution in the [Oil] and H5 features from low-to-moderate redshift. However, there is no 
obvious trend to higher redshifts. At z ~ 0.4 — 0.9, the structures all occupy regions of phase 
space which corresponds to activity levels higher than local clusters. In all cases, they lie 
above the normal star-formation lines but with large cluster-to-cluster variations, indicating 
different evolutionary histories and, subsequently, different starburst contributions. 
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As a caveat, we note that an exact comparison of CI 0023 and RX J1821 to other 
surveys depends on the spectroscopic selection criteria. We have shown above that we can 
directly compare the populations of CI 0023 and RX J 1821. However, the DS and Morphs 
spectroscopic samples are derived largely (although not entirely; see Dressier et al. 1999) 
from magnitude-limited surveys. Because red (presumably less active) galaxies are given a 
higher priority in our spectroscopic survey, the resulting composite spectra of CI 0023 and RX 
J1821 may be biased, relative to the DS and Morphs samples, toward smaller [Oil] and H<5 
EWs. As such, our measurements would represent lower limits on the star formation activity. 
In the case of CI 0023, correcting this bias would likely increase its already substantial offset 
from the lower-redshift samples and the normal star-formation lines, further confirming the 
extreme level of activity in this system. 

To quantify this bias, we plan to measure the average (and range) of spectral properties 
(e.g., EWs) as a function of color and magnitude once the spectroscopic data are completed 
for a significant portion of our survey. Incompleteness, especially in the faint blue population, 
can be measured by statistically subtracting the typical field contribution in any given color 
and magnitude range using deep imaging taken with the same filters as part of the CFHTLS. 
Combining the two, we can estimate (and account for) the effect of bias on the cluster 
properties, in particular the blue fractions and composite spectra. In addition, we will 
use photometric redshifts where possible to more accurately determine our spectroscopic 
selection function. However, accurate photometric redshifts require mid-infrared data in the 
IRAC bands which are currently only available for the CI 1604 supercluster field. 



5. Conclusions 

In this paper, we have presented the motivation, design, and implementation of the Ob- 
servations of Redshift Evolution in Large Scale Environments (ORELSE) survey, a systematic 
search for structure on scales of 10 Mpc around 20 well-known clusters at 0.6 < z < 1.3. 
The survey covers 5 square degrees, all targeted at intermediate-to-high-density regions, 
making it complementary and comparable to field surveys such as DEEP2 and COSMOS. 
The program utilizes optical/near- infrared imaging at the UKIRT 3.8-m, KPNO 4-m, Palo- 
mar 5-m, and Subaru 8-m covering at least 25' around each target cluster. Following the 
successful application in the CI 1604 supercluster at z ~ 0.9, we use an adaptively smoothed 
red-galaxy density map to visually identify associated groups/clusters and larger scale fila- 
ments/walls. Guided largely by actual observations, we use the extensive sample of 400+ 
confirmed members in CI 1604 supercluster to adapt our color and magnitude cuts, quantify 
significant density peaks, and estimate false contamination rates for other target fields at 
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lower and higher redshifts. We find that this technique is highly efficient at detecting sys- 
tems even down to group masses (a ~ 300 km s -1 ), as well as extended structures covering 
significant portions of the imaging field. 

The crucial component of the ORELSE survey, and what distinguishes it from similar 
surveys, is the unprecedented amount of spectroscopic coverage. Utilizing the wide field, 
multi-object spectrograph DEIMOS on the Keck 10-m, we are obtaining high quality spectra 
for 100-200+ confirmed members per system, allowing us to measure properties on a galaxy- 
to-galaxy basis. Targeting galaxies down to %' = 24, our system of color selection provides a 
spectroscopic efficiency for cluster members of up to 45%. 

Based on the first results from the survey, we already see significant diversity of large- 
scale and galaxy-scale properties at a given redshift. Although at similar redshifts, CI 0023 
and RX J1821 differ in their evolutionary state, global dynamics, activity level, and galaxy 
population. The optically-selected CI 0023 shows multiple high-density peaks in the red- 
galaxy density map and a complex velocity structure, suggesting a four-way group-group 
merger. The measured velocity dispersions of the member groups range from 206-479 km 
s -1 . The overall galaxy population is exceptionally active, with a high fraction of blue 
galaxies and [Oil] emitters. Conversely, the X-ray selected RX J1821 is a relatively isolated, 
centrally-condensed massive cluster. Only one significant overdensity and one redshift peak 
are observed in this field. The measured velocity dispersion is 926 km s _1 , indicative of a rich 
cluster. However, there are obvious signs of continuing dynamical evolution, including an 
elongated galaxy distribution and significant velocity substructure. The galaxy population 
is substantially more quiescent than that of CI 0023; however, it still has a larger fraction 
of blue, star-forming galaxies than local clusters, consistent with the bluing of the cluster 
population with redshift. Composite spectra from both systems suggest some contribution 
from galaxies which have undergone a recent starburst (strong H<5 absorption), with the 
contribution in CI 0023 being substantial. 

These two systems indicate that global cluster properties (i.e., mass) play a role in 
driving galaxy evolution (see also Discussion of Poggianti et al. 2006). CI 0023 is comprised 
of group-sized halos, has a small quiescent population, and subsequently contains more active 
galaxies. CI 0023 will likely be a massive cluster, with a velocity dispersion on the order of 
RX J1821, in ~ 1 Gyr, although noticeable dynamical evolution will still be taking place 
after ~ 3 Gyr, i.e. until z < 0.3 (Lubin, Postman & Oke 1998). If this dynamical evolution 
does take place, the star formation in a significant fraction of the galaxy population must be 
quenched during this period in order to be consistent with average star-formation rates and 
blue fractions in moderate- redshift clusters of similar mass (e.g., Balogh et al. 1997; Nakata et 
al. 2005). The current activity level in CI 0023 and this evolutionary time frame implies that 
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the group and/or merger environment is efficient at inducing and quenching star formation. 
Conversely, RX J1821 is already a cluster-sized halo, with a galaxy population dominated by 
passive red galaxies (from both early-epoch formation and later-time quenching). Although 
both systems will end up as massive clusters by the present day, they have clearly different 
evolutionary histories. We note that these results, along with preliminary findings from 
other ORELSE targets, suggest that X-ray-selection may be biased, at least with respect 
to optical surveys, toward more evolved systems at a given redshift. On the other hand, 
due to projection effects, optical surveys detect multi-group systems and systems embedded 
in large scale structures that are not traditionally part of "cluster" surveys. As a result, 
the cluster detection technique may have a strong impact on any conclusions concerning the 
presence and effect of large scale structure, the relative importance of global versus local 
environment, and the timescales for galaxy and cluster evolution. 

The scientific goals of the ORELSE survey are to identify and examine a statistical 
sample of dynamically active clusters and large scale structures during an active period in 
their history. Based on our completed (including spectroscopic) observations so far, such as 
those presented here and in Gal et al. (2008), we have found two superclusters (containing 
7+ groups/clusters), a four- way group-group merger, and two largely-isolated, massive X- 
ray-selected clusters. When the survey is complete, we will have a significant sample with 
which to constrain large scale cluster dynamics, cluster formation mechanisms, and the effect 
of environment on galaxy evolution. Combining our ground-based databases with follow-up 
radio, X-ray, and/or mid-infrared observations, the program will (1) determine the relation 
between global cluster properties (such as richness, mass, luminosity, and dynamical state) 
and the nature of the surrounding large scale structures; (2) quantify the dependence of 
galaxy multi-wavelength and spectral properties on position and density to determine when, 
where, and why star formation and nuclear activity are triggered (or truncated) in and near 
the cluster environment; and (3) chart the process of structure formation and its effect on 
galaxy evolution over this active 3 Gyr period. 
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manuscript. This material is based upon work supported by the National Aeronautics and 
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the W.M. Keck Foundation. The authors wish to recognize and acknowledge the very signif- 
icant cultural role and reverence that the summit of Mauna Kea has always had within the 



-26- 



indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct 
observations from this mountain. The observing staff, telescope operators, and instrument 
scientists at Keck provided a great deal of assistance. 

REFERENCES 

Andreon, S., Valtchanov, L, Jones, L. R., Altieri, B., Bremer, M., Willis, J., Pierre, M., & 
Quintana, H. 2005, MNRAS, 359, 1250 

Balogh, M. L., Morris, S. L., Yee, H. K. C, Carlberg, R. G., & Ellingson, E. 1997, ApJ, 488, 
L75 

Balogh, M. L., Morris, S. L., Yee, H. K. C, Carlberg, R. G., & Ellingson, E. 1999, ApJ, 527, 
54 

Balogh, M., et al. 2004, MNRAS, 348, 1355 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Best, P. N., van Dokkum, P. G., Franx, M., Rottgering, H. J. A. 2002, MNRAS, 330, 17 
Best, P. N. 2002, MNRAS, 336, 1293 
Best, P. 2003, Ap&SS, 285, 137 
Blakeslee, J. P., et al. 2006, ApJ, 644, 30 
Blanton, M. R., et al. 2003, ApJ, 594, 186 

Blanton, E. L., Gregg, M. D., Helfand, D. J., Becker, R. H., & White, R. L. 2003, AJ, 125, 
1635 

Bower, R. G., & Balogh, M. L. 2004, Clusters of Galaxies: Probes of Cosmological Structure 
and Galaxy Evolution, 325 

Branchesi, M., Gioia, I. M., Fanti, C, Fanti, R., & Perley, R. 2006, A&A, 446, 97 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Butcher, H., & Oemler, A., Jr. 1984, ApJ, 285, 426 

Cappelluti, N., Cappi, M., Dadina, M., Malaguti, G., Branchesi, M., D'Elia, V., & Palumbo, 
G. G. C. 2005, A&A, 430, 39 



Chang, R., Gallazzi, A., Kauffmann, G., Chariot, S., Ivezic, Z., Brinchmann, J., & Heckman, 
T. M. 2006, MNRAS, 366, 717 

Cohn, J. D., & White, M. 2005, Astroparticle Physics, 24, 316 

Colberg, J. M., et al. 2000, MNRAS, 319, 209 

Colless, M., et al. 2001, MNRAS, 328, 1039 

Davis, M., Tonry, J., Huchra, J., & Latham, D. W. 1980, ApJ, 238, L113 
Davis, M., et al. 2003, Proc. SPIE, 4834, 161 
De Lucia, C, et al. 2004, ApJ, 610, L77 
De Lucia, C, et al. 2007, MNRAS, 374, 809 
Demarco, R., et al. 2007, ApJ, 663, 164 

Dolag, K., Meneghetti, M., Moscardini, L., Rasia, E., & Bonaldi, A. 2006, MNRAS, 370, 
656 

Dressier, A., & Shectman, S. A. 1988, A J, 95, 985 
Dressier, A., et al. 1997, ApJ, 490, 577 

Dressier, A., Smail, I., Poggianti, B. M., Butcher, H., Couch, W. J., Ellis, R. S., & Oemler, 
A. J. 1999, ApJS, 122, 51 

Dressier, A., Oemler, A. J., Poggianti, B. M., Smail, I., Trager, S., Shectman, S. A., Couch, 
W. J., & Ellis, R. S. 2004, ApJ, 617, 867 

Ebeling, H., Barrett, E., & Donovan, D. 2004, ApJ, 609, L49 

Eastman, J., Martini, P., Sivakoff, C, Kelson, D. D., Mulchaey, J. S., & Tran, K.-V. 2007, 
ApJ, 664, L9 

Eke, V. R., Navarro, J. F., & Frenk, C. S. 1998, ApJ, 503, 569 

Einasto, M., Tago, E., Jaaniste, J., Einasto, J., & Andernach, H. 1997, A&AS, 123, 119 
Ellingson, E., Lin, H., Yee, H. K. C, & Carlberg, R. G. 2001, ApJ, 547, 609 
Evrard, A. E. et al. 2002, ApJ, 573, 7 
Faber, S. M. et al. 2003, Proc. SPIE, 4841, 1657 



-28- 

Fadda, D., Biviano, A., Marleau, F. R., Storrie-Lombardi, L. J., & Durret, F. 2008, ApJ, 
672, L9 

Fassbender, R. et al. 2008, A&A, in press 

Finn, R. A., Zaritsky, D., & McCarthy, D. W., Jr. 2004, ApJ, 604, 141 
Fisher, D., Fabricant, D., Franx, M., & van Dokkum, P. 1998, ApJ, 498, 195 
Frenk, C. S., Evrard, A. E., White, S. D. M., & Summers, F. J. 1996, ApJ, 472, 460 
Gal, R. R. & Lubin, L. M. 2004, ApJ, 607, LI 

Gal, R. R., de Carvalho, R. R., Lopes, P. A. A., Djorgovski, S. G., Brunner, R. J., Mahabal, 
A., & Odewahn, S. C. 2003, AJ, 125, 2064 

Gal, R. R., Lubin, L. M., & Squires, G. K. 2005, AJ, 129, 1827 

Gal, R. R., Lemaux, B. C, Lubin, L. M., Kocevski, D. & Squires, G. K. 2008, ApJin press 
Geller, M. J., & Huchra, J. P. 1989, Science, 246, 897 

Gilbank, D. G., Yee, H. K. C, Ellingson, E., Hicks, A. K., Gladders, M. D., Barrientos, 
L. F., & Keeney, B. 2008, ApJ, 677, L89 

Gilmour, R., Gray, M. E., Almaini, O., Best, P., Wolf, C, Meisenheimer, K., Papovich, C, 
& Bell, E. 2007, MNRAS, 380, 1467 

Gioia, I. M., Henry, J. P., Mullis, C. R., Ebeling, H., & Wolter, A. 1999, AJ, 117, 2608 

Gioia, I. M., Henry, J. P., Mullis, C. R., Bohringer, H., Briel, U. G., Voges, W., & Huchra, 
J. P. 2003, ApJS, 149, 29 

Gioia, I. M., Wolter, A., Mullis, C. R., Henry, J. P., Bdhringer, H., & Briel, U. G. 2004, 
A&A, 428, 867 

Girardi, M., Giuricin, G., Mardirossian, F., Mezzetti, M., & Boschin, W. 1998, ApJ, 505, 74 

Gladders, M. D., & Yee, H. K. C. 2000, AJ, 120, 2148 

Gladders, M. D., & Yee, H. K. C. 2005, ApJS, 157, 1 

Gladders, M. D., Yee, H. K. C, & Ellingson, E. 2002, AJ, 123, 1 

Gomez, P. L., et al. 2003, ApJ, 584, 210 



-29- 

Goto, T., et al. 2003a, PASJ, 55, 757 

Goto, T., Yamauchi, C, Fujita, Y., Okamura, S., Sekiguchi, M., Smail, I., Bernardi, M., & 
Gomez, P. L. 2003b, MNRAS, 346, 601 

Gunn, J. E., & Gott, J. R. I. 1972, ApJ, 176, 1 

Gunn, J. E., Hoessel, J. G., & Oke, J. B. 1986, ApJ, 306, 30 

Gray, M. E., Taylor, A. N., Meisenheimer, K., Dye, S., Wolf, C, & Thommes, E. 2002, ApJ, 
568, 141 

Hashimoto, Y., Henry, J. P., Hasinger, G., Szokoly, G., & Schmidt, M. 2005, A&A, 439, 29 

Henry, J. P., Mullis, C. R., Voges, W., Bohringer, H., Briel, U. G., Gioia, I. JVL, & Huchra, 
J. P. 2006, ApJS, 162, 304 

Holden, B. P., Stanford, S. A., Squires, G. K., Rosati, P., Tozzi, P., Eisenhardt, P., & 
Spinrad, H. 2002, AJ, 124, 33 

Homeier, N. L., et al. 2006, ApJ, 647, 256 

Howell, S. B., & Tonry, J. L. 2003, Bulletin of the American Astronomical Society, 35, 1235 

Jannuzi, B. T., et al. 2004, Bulletin of the American Astronomical Society, 36, 1478 

Jannuzi, B. T., & Dey, A. 1999, Photometric Redshifts and the Detection of High Redshift 
Galaxies, 191, 111 

Jeltema, T. E., Canizares, C. R., Bautz, M. W., & Buote, D. A. 2005, ApJ, 624, 606 

Kaiser, N., Wilson, G., Luppino, G., & Dahle, H. 1999, ArXiv Astrophysics e-prints, 
a?Xiv:astro-ph/9907229 

Kartaltepe, J.S., Ebeling, H., Ma, C.J. & Donovon, D. 2008, MNRAS, submitted 

Kocevski, D., et al. 2008a, ApJ, in press 

Kocevski, D., et al. 2008b, ApJ, submitted 

Kodama, T., Smail, I., Nakata, F., Okamura, S., & Bower, R. G. 2001, ApJ, 562, L9 
Kodama, T., Smail, I., Nakata, F., Okamura, S., & Bower, R. G. 2003, ApJ, 591, L169 
Kodama, T., et al. 2005, PASJ, 57, 309 



-30- 

Kodama, T., & Arimoto, N. 1997, A&A, 320, 41 
Koekemoer, A. M., et al. 2007, ApJS, 172, 196 
Koester, B. P., et al. 2007, ApJ, 660, 221 

Koyama, Y., Kodama, T., Tanaka, M., Shimasaku, K., & Okamura, S. 2007, MNRAS, 382, 
1719 

Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 
Lewis, L, et al. 2002, MNRAS, 334, 673 

Lopes, P. A. A., de Carvalho, R. R., Capelato, H. V., Gal, R. R., Djorgovski, S. G., Brunner, 
R. J., Odewahn, S. C, & Mahabal, A. A. 2006, ApJ, 648, 209 

Lubin, L. M., Postman, M., Oke, J. B., Ratnatunga, K. U., Gunn, J. E., Hoessel, J. G., & 
Schneider, D. P. 1998, AJ, 116, 584 

Lubin, L. M., Postman, M., & Oke, J. B. 1998, AJ, 116, 643 

Lubin, L. M., Brunner, R., Metzger, M. R., Postman, M., & Oke, J. B. 2000, ApJ, 531, L5 

Lubin, L. M., Oke, J. B., & Postman, M. 2002, AJ, 124, 1905 

Lubin, L. M., Mulchaey, J. S., & Postman, M. 2004, ApJ, 601, L9 

Margoniner, V. E., & Wittman, D. M. 2008, ApJ, 679, 31 

Marcillac, D., Rigby, J. R., Rieke, G. H., & Kelly, D. M. 2007, ApJ, 654, 825 

Marcillac, D., et al. 2008, ApJ, 675, 1156 

Maughan, B. J., Jones, L. R., Ebeling, H., & Scharf, C. 2006, MNRAS, 365, 509 
Mei, S., et al. 2006a, ApJ, 644, 759 
Mei, S., et al. 2006b, ApJ, 639, 81 
Mihos, J. C. 1995, ApJ, 438, L75 
Mihos, C. 1999, Ap&SS, 266, 195 

Mihos, J. C. 2004, Clusters of Galaxies: Probes of Cosmological Structure and Galaxy 
Evolution, 277 



-31 - 

Miller, N. A., & Owen, F. N. 2003, AJ, 125, 2427 

Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 

Mullis, C. R., et al. 2003, ApJ, 594, 154 

Nakata, F., et al. 2005, MNRAS, 357, 1357 

Oke, J. B., Postman, M., & Lubin, L. M. 1998, AJ, 116, 549 

Oke, J. B., et al. 1995, PASP, 107, 375 

Owen, F. N., Ledlow, M. J., Keel, W. C, & Morrison, G. E. 1999, AJ, 118, 633 

Patel, S. G., Kelson, D. D., Holden, B. P., Illingworth, G. D., Franx, M., van der Wei, A., 
& Ford, H. 2008, iarXiv:0812.202T1 

Pierre, M., et al. 2006, MNRAS, 372, 591 

Poggianti, B. M., et al. 2006, ApJ, 642, 188 

Porter, S. C, Raychaudhury, S., Pimbblet, K. A., & Drinkwater, M. J. 2008, MNRAS, 388, 
1152 

Postman, M., Geller, M. J., & Huchra, J. P. 1988, AJ, 95, 267 
Postman, M., Lubin, L. M., & Oke, J. B. 1998, AJ, 116, 560 

Postman, M., Lubin, L. M., Gunn, J. E., Oke, J. B., Hoessel, J. G., Schneider, D. P., & 
Christensen, J. A. 1996, AJ, 111, 615 

Quintana, H., Ramirez, A., Melnick, J., Raychaudhury, S., & Slezak, E. 1995, AJ, 110, 463 

Rines, K., Geller, M. J., Diaferio, A., Mahdavi, A., Mohr, J. J., & Wegner, G. 2002, AJ, 
124, 1266 

Scoville, N., et al. 2007, ApJS, 172, 1 

Simcoe, R. A., Metzger, M. R., Small, T. A., & Araya, G. 2000, Bulletin of the American 
Astronomical Society, 32, 758 

Smail, I., Edge, A. C, Ellis, R. S., & Blandford, R. D. 1998, MNRAS, 293, 124 

Small, T. A., Ma, C.-P, Sargent, W. L. W., & Hamilton, D. 1998, ApJ, 492, 45 

Stanford, S. A., Eisenhardt, P. R., & Dickinson, M. 1998, ApJ, 492, 461 



-32 - 



Stanford, S. A., Holden, B., Rosati, P., Eisenhardt, P. R., Stern, D., Squires, G., & Spinrad, 
H. 2002, AJ, 123, 619 

Swinbank, A. M., et al. 2007, MNRAS, 379, 1343 

Tanaka, M., Goto, T., Okamura, S., Shimasaku, K., & Brinkmann, J. 2004, AJ, 128, 2677 

Tanaka, M., Kodama, T., Arimoto, N., Okamura, S., Umetsu, K., Shimasaku, K., Tanaka, 
L, & Yamada, T. 2005, MNRAS, 362, 268 

Tanaka, M., Kodama, T., Arimoto, N., & Tanaka, I. 2006, MNRAS, 365, 1392 

Tanaka, M., Hoshi, T., Kodama, T., & Kashikawa, N. 2007a, MNRAS, 379, 1546 

Tanaka, M., Kodama, T., Kajisawa, M., Bower, R., Demarco, R., Finoguenov, A., Lidman, 
C, & Rosati, P. 2007b, MNRAS, 377, 1206 

Tody, D. 1986, Proc. SPIE, 627, 733 

Tran, K.-V. H., Franx, M., Illingworth, G., Kelson, D. D., & van Dokkum, P. 2003, ApJ, 
599, 865 

Treu, T., Ellis, R. S., Kneib, J. -P., Dressier, A., Smail, I., Czoske, O., Oemler, A., & 
Natarajan, P. 2003, ApJ, 591, 53 

Valtchanov, I., et al. 2004, A&A, 423, 75 

van Dokkum, P. G., & Franx, M. 2001, ApJ, 553, 90 

van Dokkum, P. G., Franx, M., Fabricant, D., Illingworth, G. D., k Kelson, D. D. 2000, 
ApJ, 541, 95 

van Dokkum, P. G., Stanford, S. A., Holden, B. P., Eisenhardt, P. R., Dickinson, M., & 
Elston, R. 2001, ApJ, 552, L101 

Vikhlinin, A., McNamara, B. R., Forman, W., Jones, C, Quintana, H., & Hornstrup, A. 
1998, ApJ, 498, L21 

Yan, R., Newman, J. A., Faber, S. M., Konidaris, N., Koo, D., & Davis, M. 2006, ApJ, 648, 
281 

Zabludoff, A. I., & Mulchaey, J. S. 1998, ApJ, 496, 39 



This preprint was prepared with the AAS IATgX macros v5.2. 



Table 1. The Cluster Sample 



Cluster 


R.A. (J2000) 


Dec. (J2000) 


2 


Selection 


Chandra/XMM 
Observations? 


Optical 
Complete? 


Near-infrared 
Complete? 


Notes 


Refs 


pi nnoQj_n/ioQ 


00 23 


52.2 


+04 


23 


07 


U.O40 


Optical 


1 Ob 


Vac 
1 Ob 


Partial 


complex merger 


(1J 


r pq Tnoo/i nnno 


02 24 


30.0 


—00 


02 


00 


n 770 
V. 1 I z 


Optical 


1 es 






strong lens 


+ J 




02 27 


09 7 


—04 


18 


05 


1 nnn 
l.UUU 


X-ray 


Vac 
I Ob 


Yes 




double cluster 


(■>) 


Kji JUo4y-t-440Z 


08 48 


56 3 


+44 


52 


16 


1 .ZOl 


X-ray 


1 Ob 


1 es 




Lynx double cluster 




PI TnQ1fl-i-^499 


09 10 


44.9 


-1-54 


_»j 


09 


1 1 1 n 


X-ray 


1 es 


Voo 

1 es 


1 (. h 


elongated structure 




PI flQQ/1 Sfl/1 


09 43 


41.0 


+48 


04 


46 


u.oyy 


Optical 


1 Ob 


Vac 
1 Ob 


1 es 






xl-A- J 1UDO. 1 to ( oO 


10 53 


39.8 


+57 


35 


18 


1 1 a n 


X-ray 


Vac 

1 Ob 


i es 


Vac 

i eb 


double cluster 


(>) 


llo 1 tOUUU 


11 37 


33.4 


+30 


07 


36 


u.yoy 


Radio 




1 es 


1 Cb 


wide-angle-tailed radio source 




RX J1221. 4+4918 


12 21 


24.5 


+49 


18 


13 


0.700 


X-ray 


Yes 


Yes 


Yes 


double cluster 


(9) 


CI 1324+3011 


13 24 


50.4 


+30 


11 


26 


0.756 


Optical 


Yes 


Yes 


Yes 


supercluster member 


(10) 


CI 1324+3059 


13 24 


53.5 


+30 


59 


00 


0.755 


Optical 


Yes 


Yes 


Yes 


supcrcluster member 


(11) 


CI 1325+3009 


13 25 


18.7 


+30 


09 


57 


1.067 


Optical 




Yes b 


Yes 




(12) 


CI J1350+6007 


13 50 


48.5 


+60 


07 


07 


0.804 


X-ray 


Yes 


Yes 


Yes 


disturbed, elliptical structure 


(13) 


CI J1429. 0+4241 


14 29 


06.4 


+42 


41 


10 


0.920 


X-ray 


Yes 


Yes 


Yes 


elongated structure 


(14) 


CI 1604+4304 


16 04 


19.5 


+43 


04 


33 


0.900 


Optical 


Yes 


Yes 


Yes 


supcrcluster member 


(15) 


CI 1604+4321 


16 04 


31.5 


+43 


21 


17 


0.920 


Optical 


Yes 


Yes 


Yes 


supercluster member 


(16) 


RX J1716. 4+6708 


17 16 


49.6 


+67 


08 


30 


0.813 


X-ray 


Yes 


Yes b 


Yes 


elongated core, small subclustcr 


(17) 


RX J1757.3+6631 


17 57 


19.4 


+66 


31 


31 


0.691 


X-ray 


Yes 


Yes 


Yes 




(18) 


RX J1745. 2+6556 


17 45 


18.2 


+65 


55 


42 


0.608 


X-ray 




Yes 


Yes 




(19) 


RX J1821.6+6827 


18 21 


32.9 


+68 


27 


55 


0.818 


X-ray 


Yes 


Yes 


Yes 


elongated structure 


(20) 



CO 
CO 



a Arhival data from the CFHT Legacy Survey (see http://www.cfht.hawaii.edu/Sclence/CFHTLS/ ). 
b Arhival data from Subaru (see http://smoka.nao.ac.jp/ ). 
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& Oke (1986); Oke, Postman, & Lubin (1998); (12) Postman et al. (1996); Oke, Postman, & Lubin (1998); (13) Rosati et al. (1998); Holden et al. (2002); 
Jeltema et al. (2005); (14) Maughan et al. (2006); (15) Gunn, Hoessel, & Oke (1986); Oke, Postman, & Lubin (1998); (16) Gunn, Hoessel, & Oke (1986); Oke, 
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Table 2. The Red-Galaxy Color Selection 



Rcdshift Range 


r' - 


-%' 


i' - 


z' 


0.6 - 0.7 


1.0 - 


- 1.4 


0.3 - 


0.7 


0.7 - 0.8 


1.0 - 


- 1.4 


0.4- 


0.8 


0.8 - 0.9 


1.0 - 


- 1.4 


0.5 - 


0.9 


0.9 - 1.0 


1.0 - 


- 1.4 


0.6 - 


1.0 


> 1.0 






> 0.9 
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Table 3. Target Coordinates, Redshifts, and Velocity Dispersions 



J2000.0 Within 0.5h^g Mpc Within l.Ohfp 1 Mpc 



Name 


ID 


RA 


Dec 


N 


Zmean 




(7 




N 










CI 0023+0423 


A 


6.0223 


4.3586 


5 


0.8383 


a 






7 


0.8390 


479 


± 


172 




Bl b 


5.9764 


4.3874 


10 


0.8286 


242 


± 


54 


17 


0.8282 


206 


± 


32 




B2 b 


5.9671 


4.3824 


12 


0.8445 


177 


± 


118 


27 


0.8453 


293 


± 


79 




C 


5.9227 


4.3824 


15 


0.8451 


548 


± 


84 


29 


0.8463 


428 


± 


67 


RX J1821.6+6827 




275.3801 


68.4651 


17 


0.8170 


993 


± 


87 


10 


0.8193 


926 


± 


77 



insufficient redshifts to compute dispersion 

b There are two distinct velocity peaks in this component (see Figure IS}. 
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Fig. 1. — The r' — i' vs. i' — z' color-color plane. The evolution with redshift for three 
solar-metallicity BC03 models are shown (see §2.3p as dotted lines with the colors at a given 
redshift indicated. Large filled squares represent the average colors of spectroscopically- 
confirmed red-sequence galaxies in three ORELSE structures, the CI 1604 supercluster at 
z m 0.9 (cyan), the CI 1324 supercluster at z w 0.73 (magenta), and RX J1821 at z = 0.82 
(green). Solid lines show the offset from the observed colors of these three systems to the 
colors predicted by the instantaneous burst model at Zf = 3. No single model fits all three 
systems. 
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2 




z' r'-i' 

Fig. 2. — The i' vs. (r' — i'), z' vs. (r' — z'), z' vs. (i' — z'), and (r' — i') vs. (z' — z') color- 
magnitude and color-color diagrams of the CI 0023 field. All objects in the LFC imaging are 
shown as small dots. The black rectangular regions outline the color-selections applied to 
produce the density maps and prioritize spectroscopic targets. Spectroscopically confirmed 
system members are overplotted as large red dots. Foreground galaxies at z < 0.82 are 
shown as blue dots, background galaxies at z > 0.87 are green dots, and cyan dots are stars. 
Color distributions of spectroscopic objects are shown to the right of each CMD. 
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Fig. 3. — Same as Figure [2] but for the RX J1821 field. Spectroscopically confirmed cluster 
members are overplotted as large red dots. Foreground galaxies at z < 0.80 are shown as 
blue dots, background galaxies at z > 0.84 are green dots, and cyan dots are stars. 
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Fig. 4. — Adaptive kernel density map of color-selected red galaxies in the CI 0023 field. 
(Left) Candidate groups/clusters, which meet our detection threshold (see §2.41) . are marked 
with red circles of 0.5 Mpc and labeled with an identifying letter. Large blue rectangles 
indicate the layout of the five DEIMOS masks. (Right) Small circles indicate the 169 con- 
firmed members, coded by color representing four redshift ranges which correspond to peaks 
in the redshift distribution in Figure [61 The distribution of confirmed members corresponds 
extremely well to peaks in the density map, including, as expected, the three significant 
peaks but also the lower-density peak directly South of component B. 
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Fig. 5. — Adaptive kernel density map of color-selected red galaxies in the RX J1821 field. 
(Left) Only the cluster itself, marked by a red circle of radius 0.5 Mpc, meets our 
detection threshold. Large blue rectangles indicate the layout of the two DEIMOS masks. 
(Right) Small red circles indicate the 85 confirmed cluster members. As in Figure HI the 
distribution of confirmed members follows extremely well the overdense regions in the density 
map, tracing the elongated structure of the cluster and even picking up the lower-density 
peaks (small groups) to the South and South-west. 
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Fig. 6. — Redshift distribution in the CI 0023 field. The top panel shows all extragalactic 
objects (dotted line) and only those with high redshift quality (solid line). The bottom panel 
focuses on the redshift range of the group-group merger with redshift bins of Az = 0.001. 
There are four clear peaks in the redshift histogram, indicative of the dynamical state of the 
system. 
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Fig. 7. — Same as Figure [6] for the RX J1821 field. There is a single clear peak z ~ 
0.817 which is completely isolated in redshift space; however, the distribution shows clear 
indications of velocity substructure (see also Figure [9]). 
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Fig. 8. — Redshift histograms for components A, B, and C identified in the density map of 
the CI 0023 field (see Figure 0J. The shaded regions show the distribution within a projected 
radius of 0.5 h^Q Mpc, while the solid line corresponds to a radius of 1 h^l Mpc. Components 
A and C have reasonably well-defined redshift peaks with measured velocity dispersions that 
are consistent with galaxy groups (see TableE]), while component B is actually a superposition 
of two distinct groups along the line of sight. 
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Fig. 9. — Velocity histograms of the cluster RX J1821 within projected radii of 0.5, 1.0, 
and 1.5 h^Q Mpc. Each velocity is plotted relative to the median velocity of the member 
galaxies in each bin. Although the velocity dispersion estimators are quite consistent within 
a radial bin and from bin to bin, the distributions are clearly non-gaussian. The cluster does 
not, however, have clear multiple components like CI 0023 or a significant amount of nearby 
structure like the CI 1604 supercluster. 
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Fig. 10. — Histogram of r' — i' (left panel) and i' — z' (right panel) colors of confirmed member 
galaxies (down to %' = 24) in CI 0023 (dotted histogram) and RX J1821 (solid histogram). 
RX J 1821 shows a tighter red-sequence and a smaller percentage of blue galaxies than CI 
0023. 
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Fig. 11. — A comparison between the composite DEIMOS spectra of CI 0023 (solid black 
line) and RX J1821 (dashed magenta line). Only confirmed members with redshift quality 
Q > 3 are used in both cases: 138 and 69 members for CI 0023 and RX J1821, respectively 
(see §4.21) . The two composites are normalized so that their fluxes match at 4050 A. Insets 
in the lower right show close ups of the [Oil] and B.5 lines. As already indicated by their 
colors, the galaxy population of CI 0023 is considerably more active than that of RX J1821. 
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Fig. 12. — Adapted Dressier et al. (2004) plot of [Oil] vs. H5 equivalent width. The two 
dashed lines (based on the CNOC2 and 2dF surveys) indicate the average spectral properties 
of a galaxy population consisting of passive galaxies plus a varying fraction of continuously 
star-forming galaxies. The small asterisks along the curves mark the appropriate values 
of EWQOII]) and EW(FW) for 20%, 40%, 60%, 80%, and 100% mix of galaxies undergoing 
continuous star formation. Local clusters fall on the lines, while mo derate- redshift clusters, as 
does the massive cluster MS 1054 at z — 0.83 (open green square), lie above indicating some 
contribution from galaxies that have undergone a recent starburst. The average properties of 
the CI 0023 system (large black circle) is substantially offset in both [Oil] and H<5, implying a 
high fraction of star-forming galaxies and an even larger contribution from recent starbursts. 
In contrast, the X-ray-selected, isolated cluster RX J1821 at z = 0.82 (large black diamond) 
has more moderate [Oil] emission and falls close to the normal star-formation lines, indicating 
a large quiescent galaxy population (as also evident in the CMD of Figure [3]) and only a 
modest starburst contribution. 



